Astronomy & Astrophysics manuscript no. belloche 


February 2, 2008 


(DOI: will be inserted by hand later) 





Disappearance of N 2 H + from the Gas Phase 
in the Class Protostar IRAM 04191* 



A. Belloche 12 and P. Andre 3 



o 
o 

(N 

Sh 

Oh 
< 

in 



> 

Os 
O 

o 

Oh 
i 

o 

c3 



1 Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
e-mail: belloche@mpifr-bonn.mpg.de 

2 LERMA/LRA, Ecole Normale Superieure, 24 rue Lhomond, F-75231 Paris Cedex 05, France 

3 Service d'Astrophysique, CEA/DSM/DAPNIA, C.E. Saclay, F-91191, Gif-sur-Yvette Cedex, France 
e-mail: pandre@cea.fr 

Received 23 February 2004; accepted 10 April 2004 

Abstract. We present a high-resolution millimeter study of the very young Class protostar IRAM 04191 + 1522 in the Taurus 
molecular cloud. N2H + (l-0) observations with the IRAM Plateau de Bure Interferometer and 30m telescope demonstrate that 
the molecular ion N 2 H + disappears from the gas phase in the inner part of the protostellar envelope (r < 1600 AU, tih, > 5 x 10 5 
cur 3 ). This result departs from the predictions of current chemical models. It suggests either that N 2 is more depleted than the 
models predict, owing to a higher binding energy on polar ice or an enhanced grain chemistry transforming N 2 to less volatile 
species, or that strong deuterium fractionation enhances N 2 D + to the detriment of N 2 H + . 
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1. Introduction 

Understanding the onset of gravitational collapse in dense 
cloud cores requires detailed (sub)millimeter studies of the 
structure of preste llar condensations and young protostars (e.g. 
lAndre et"afll200 0). However, recent observations have shown 



that molecules such as CO and CS deplete onto grai n sur- 
faces in the inner parts of dense cores (e.g. iBacmann et alJ 
l2002HTafalla et alJl2002h . Since this depletion phenomenon is 
thought to affect many other species at high densities, study- 
ing the kinematics of pre/protostellar cores with future high- 
resolution instruments such as ALMA may be difficult and re- 
quires the identification of the best tracers of dense gas. N2H + 
has been put forward as such a good tracer since it is much less 
sensitive to depletion effects tha n other species observed with 
single-dish telescopes (e.g. Tafal la et alJl2002b . even though 
Berg in et alj f2002) reported a slight decrease of its abundance 
toward the prestellar core B68. 

Prime targets for observational studies of protostel- 
lar collapse are very young Class objects such as 
IRAM 041 91 + 1522 - hereaft er IRAM 04191 for short - in 
Taurus (see lAndre etaDl 19991 - hereafter IaMB99I) . This pro- 
tostar features a prominent (~ 0.5 - 1.5 M ) envelope and a 
powerful bipolar outflow, but lacks a sizeable accretion disk. 



Send offprint requests to: A. Belloche 

* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



iBelloche et ai1(l2002l - hereafter|BADB02) showed that the en- 
velope is undergoing both extended infall motions and fast, dif- 
ferential rotation. They proposed that the rapidly rotating inner 
envelope (r < 3500 AU) corresponds to a magnetically super- 
critical core decoupling from an environment still supported by 
magnetic fields and strongly affected by magnetic braking. 

Here, we report new observations of IRAM 04191 car- 
ried out with the Plateau de Bure interferometer (PdBI) in the 
N2H + (l-0) line in an effort to probe the inner structure of the 
envelope. We discuss the results of these high-resolution obser- 
vations which show that N2H + disappears from the gas phase 
in the inner part of the envelope. 

2. Observations and data reduction 

We observed IRAM 04191 with the PdBI for 34 h between 
December 2001 and September 2002 with 4 to 6 antennas in 
the C and D configurations. In the 3mm band, each receiver 
was tuned to the N 2 H + (101-012) line at 93.176258 GHz in 
single side-band mode. The system temperatures were typi- 
cally 110-200 K. The 39 kHz channel spacin g resulted in a 
velocity resolution of 0.20 km s _1 (see rWiesemeveil20 01 ). The 
(naturally-weighted) synthesized half-power beam width was 
5.8"x4.0" (810 AU x 560 AU) and the (FWHM) primary beam 
~ 54". Several nearby phase calibrators were observed to de- 
termine the time-dependent complex antenna gains. The corre- 
lator bandpass was calibrated on the sources 3C84, 0415+379, 
3C273, 3C454.3 or NRAO150, while the absolute flux density 
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Fig. 1. N 2 H + (l-0) integrated intensity map of the IRAM 04191 
inner envelope obtained with PdBI (uncorrected for primary 
beam attenuation). The emission is integrated over three ve- 
locity intervals covering the whole seven-component multiplet. 
The contours vary from -0.5 to -0. 1 (dotted contours) and from 
0.1 to 0.6 (solid contours), by step of 0.1 Jy beam 1 km s _1 
(i.e. 3 times the rms noise level). The 54" (FWHM) primary 
beam is shown as a big circle and the synthesized clean beam 
(HPBW) is in the bottom right corner. The cross ma rks the cen - 
tral 1.3mm continuum position oflRAM 04191 (cf.lBADB02). 
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Fig. 2. N 2 H + (l-0) integrated intensity map of the IRAM 04191 
envelope combining PdBI and 30m data. The emission is inte- 
grated over the same velocity intervals as in Fig.^ The con- 
tours vary from 0.16 to 2.1 by step of 0.16 Jy beam -1 km s _1 
(i.e. 6 times the rms noise). The map is not corrected for pri- 
mary beam attenuation (automatically applied to the single- 
dish data in the combination process); the big circle represents 
the FWHM primary beam of the PdBI at 3mm. The synthesized 
clean beam (HPBW) is shown in the bottom right corner. The 
cross marks the central position of IRAM 04191. 



scale was derived from MWC349 or 3C84. The absolute cal- 
ibration uncertainty is estimated to be ~ 15%. The data were 
calibrated and imaged using the GILDAS software (Jan. 2002 
release). The deconvolution was performed with the CLEAN 
method, down to 0.6 times the rms noise level. 

The present study also uses single-dish N2H + (l-0) and 
1.3mm contin uum on-th e-fly maps taken with the I RAM 30m 
telesco pe (see iBADBoH and lMotte & AndrefeOOll - hereafter 
MA01 ), as well as 1.3mm continuum interferometric data ob- 
tained with PdBI in 1999 (see lBADB02l) . The combination of 
the PdBI and single-dish data was performed in the wv-plane 
(see|Guilloteau 2001). We used a singl e-dish calibration f actor 
of 5.9 Jy beam" 1 K" 1 for N 2 H + (l-0) JGreve et alj fl998). For 
the continuum data, we applied a calibration factor of 1.35 to 
account for the difference in observed frequency (227 GHz at 
PdBI and ~ 250 GHz for bolometer observations at the 30m). 

3. Analysis: A drop in N 2 H + abundance 

3.1. A central hole of N 2 H + emission 

The N2H + (l-0) integrated intensity map obtained with PdBI 
shows two emission peaks to the South-East and North- West 
of the central IRAM 04191 position (see Fig.Q- An emission 
gap with negative contours down to the -15cr level is also seen 
close to the central position, shifted by 3.5" (i.e., less than one 
beam width) to the North- West. These three features, which 
form a structure elongated roughly perpendicular to the outflow 
axis, have the same relative intensities in each AFi group of hy- 
perfine components of the N 2 H + (l-0) multiplet. To recover the 
extended emission filtered out by the interferometer and obtain 
reliable N 2 H + column densities, we combined our PdBI ob- 
servati ons with sh ort-spacing data obtained with the 30m tele- 
scope (BADBOI. The resulting map of N 2 H + (l-0) integrated 
intensity reveals a ring-shaped core with two local peaks, sur- 



rounding an emission gap near the center (see Fig.|2ji. This gap 
is a factor of 3.5 (2.3) weaker than the primary (secondary) 
emission peak to the South-East (West). The detection of such 
a gap strongly suggests that the N 2 H + abundance decreases to- 
ward the center. 

3.2. Optical depths and excitation temperatures 

We used the hyperfine structure fitting method of the CLASS 
software to estimate the optical depths and excitation temper- 
atures of N 2 H + (l-0) in the PdBI-30m combined map. The fits 
were limited to the spectra with a signal-to-noise ratio larger 
than 10. We fitted only the 4 (AF[ = +1) components to avoid 
convergence problems that occurred for a few positions in the 
western part of the env elope, and may be related to excitation 
anomalies (e.g. lTurnerl l2001 ) and/or low signal-to-noise ratio 
for the weakest (110-011) component. The 1-cr optical depth 
uncertainty given by the fits is typically 15-25 % for spectra 
with opacities larger than 5, and 50-150 % for optically-thinner 
spectra. Except for a few positions at ~ (0", 15"), the total 
N 2 H + (l-0) optical depths are lower than 15, and even lower 
than 10 toward the brightest emission peak and the emission 
gap, corresponding to optical depths lower than 1.7 and 1.1, 
respectively, for the isolated (101-012) component. The excita- 
tion temperature was found to vary from 4 K to 10 K across the 
PdBI-30m combined map and to display a maximum toward 
the brightest emission peak of the integrated intensity map. 

3.3. Comparison with the dust emission maps 

Figure|3]shows the ratio of N 2 H + (101-012) integrated intensity 
to 1.3mm continuum intensity across the IRAM 04191 enve- 
lope. The large gr eyscale map was obtained using the single- 
dish 1 .3mm data of MA01 smoothed to the 26.4" (FWHM) res- 
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Fig. 3. Greyscale map of the ratio of N2H + (101-012) integrated 
intensity to 1.3mm continuum intensity, both observed with 
the 30m telescope. The (thick) contours go from 6 to 18 by 
3 Jybeam -1 kms^^Jybeairr 1 ). The 1.3mm continuum map 
smoothed to the N2H + (l-0) resolution is overlaid as dashed 
contours. The insert in the top-right corner shows the PdBI- 
30m combined map of the same ratio, limited to the 22" PdBI 
primary beam at 227 GHz (corresponding to the white circle 
in the larger map). The greyscale is the same as for the larger 
map; the contours go from 3 to 9 by 3 and from 15 to 27 by 
6 Jybeam -1 kms -1 /(Jybeam -1 ). The beam sizes are shown in 
both maps. The cross marks the position of IRAM 04191, the 
star indicates the position of the Class I object IRAS 04191. 

olution of the 30m telescope at the N2H + (l-0) frequency. The 
ratio was computed only at positions where both the 1.3mm 
and N2H + (101-012) intensities were stronger than 4 times the 
rms noise. The higher-resolution, smaller greyscale map in the 
insert of Fig.|3]displays the same N2H + (101-012) to dust con- 
tinuum ratio derived from the PdBI-30m combined data for 
both the N2H + (101-012) and 1.3mm continuum emission. The 
PdBI continuum data were recalibrated to the frequency of the 
30m continuum data prior to combination (see Sect. [3, and the 
1.3mm combined map was smoothed to the resolution of the 
N2H + (l-0) combined map. The ratio of the two sets of com- 
bined data was computed after correction for the respective 
PdBI primary beam attenuations at 3mm and 1.3mm. 

The N2H + (l-0) integrated intensity is expected to primar- 
ily track the variations of N2H + column density (see Fig.^Ji). 
Indeed, N2H + (l-0) is thermalized for «h 2 > 10 5 cm -3 , and the 
gas/d ust temperatures (which are well coupled above this den- 
sity - iDotv & Neufeldlll997l) are in the range 6.5-10 K (cf. 
Fig. 7d of BADB02), i.e., in a regime where the LTE line inten- 
sities depend only weakly on temperature. The 1.3mm contin- 
uum intensity, on the other hand, should scale as the H2 column 
density times the dust te mperature times the dust opacity (see, 
e.g., Eq. 1 of M otte etall l998). The bolom etric lum inosity of 
IRAM 04191 is so low {L bM ~ 0.15 L Q - see AMB99) that the 
dust temperature is expected to vary by les s than 50% in the 
range of radii 400-2800 AU (see BADB02). The variations of 
the intensity ratio in the central r < 2800 AU (rg < 20") region 
around IRAM 04191 should thus reflect the variations of the 
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Fig. 4. a) Map of the N2H + column density A^n,h + (squares) 
derived from the N2H + (101-012) integrated intensity map 
overlaid as contours (0.03, 0.06 to 0.3 by step of 0.06 
Jy beam -1 km s -1 ). The size of each square increases linearly 
with JV N2 h+, from 2.6 x 10 12 to 3.3 x 10 13 cirT 2 . b) Circularly- 
averaged radial profile of N^h + - The crosses with error bars 
represent the average N^h* in concentric annuli limited to the 
positions shown with filled sq uares in a). The solid curve shows 
the H2 column density profile dMAOlh scaled to the 7 th annulus 
for comparison. The two dashed curves illustrate the effect of a 
central hole of N 2 H + material of radius 1350 AU and 1850 AU, 
respectively. 



N2H + abundance, to within a factor of ~ 2 if we account for 
a possible slight saturation of N2H + (101-012) (T101-012 < 1.7, 
see Sect. 13. 2> and possible changes in the dust opacity. 

The intensity ratio has an average value of 
~ 11 Jy beam -1 km s -1 /(Jybeam -1 ) in the central (rg < 20") 
region of the single-dish map, and drops to a minimum 
of 2.7 Jybeam -1 km s -1 /(Jy beam -1 ) near the center of the 
high-resolution map. We conclude that the N2FF abundance 
averaged along the line of sight decreases by a factor of ~ 4 
from r e ~ 20" to r g ~ 5". 

3.4. Disappearance of N 2 H + from the gas phase 

The map of N2FF column density shown in Fig.0t was com- 
puted using a LTE partition function at the N2H + (l-0) excita- 
tion temperature derived from the hfs fits, assumed to be uni- 
form along each line of sight. We also applied an optical depth 
correction of t/(1 - e -T ) using the hfs fit results. Figure |4j> 
compares the circularly-averaged radial profile of N2H + col- 
umn density with the H2 column density profile derived from 
the dust emission jMAOll) . The former clearly departs from the 
latter for < 13". A good model fit is obtained when one as- 
sumes a hole of N2I-U material in the central part of the en- 
velope (r < 1350-1850 AU, see dashed curves). This compar- 
ison strongly suggests that the N2FF ion disappears from the 
gas phase above a density of «h 2 ~ 4 - 7 x 10 s cm -3 in the 
IRAM 04191 envelope, according to the density profile derived 
bv lMAOll (see Fig. 7a of BADBol. We derive a mean N 2 H + 
abundance of ~ 7 ± 4 x 10" 10 at r = 3500 AU. 
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3.5. Fast, differential rotation 

We measure a large velocity difference of 0.47 + 0.05 km s _1 
between the two N2H + (l-0) emission peaks of the PdBI map 
shown in Fig^ largely exceeding the FWHM line widths mea- 
sured toward the peaks themselves (0.34+0.08 km s _1 ). A least- 
square linear fit to the LSR-velocity map of the PdBI-30m com- 
bined data yields a mean velocity gradient of ~ 17 km s _1 pc _1 
with a position angle of RA. ~ 133°. This small-scale gradi- 
ent has a direction consistent with the direction of the mean 
velocity gradient measured on larger scale with the 30m tele- 
scope by BADB02, but is 5 times larger. The PdBI results thus 
confirm and strengthen the presence of fast, differential rotation 
identified by BADB02. The detailed form of the rotation veloc- 
ity profile, whi ch is conf irmed to have a break at r i 3000 AU 
(see Fig. 12 of BADB02), will be discussed in a forthcoming 
paper based on radiative transfer simulations taking the N2FF 
depletion into account. 

4. Implications 

We showed in Sect. l3.4l that the molecular ion N2H + disappears 
from the gas phase above a density of «h, ~ 5 x 10 5 cirT 3 in the 
envelope of the Class protostar IRAM 04 1 9 1 . This conclusion 
is at variance with the results of recent ch emical models (e.g. 
IShematovich et all2003UA ikawa et all2003l) . In particular, the 
models calculated bv lAikawa et al J (120 03) including gas-phase 
reactions, gas-dust interactions, and diffusive grain-surface re- 
actions predict an increase in the N2H + abundance during the 
prestellar collapse phase up to a density of «h 2 ~ 1 x 10 7 cm -3 
or more. IRAM 04191 is so young, probably less than 3xl0 4 yr 
since the formation of the central object (AMB99), that it is dif- 
ficult to invoke a large decrease in the N2FF abundance on such 
a short timescale at the beginning of the main accretion phase. 

Since C O is one of the main destroyers of N2I-F in the gas 
phase (e.g. lAikawa et alJl200ll) . its potential desorption from 
grain mantles due to heating by the central protostar could be a 
mechanism responsible for the disappearance of N2H + from the 
gas phase. However, the luminosity of IRAM 04191 is so low 
(see Sect. 13. 3> that it is insufficient to raise the dust temperature 
above 1 1 K for r > 400 AU (J. Bouwman 2003 , private c om- 
munication). Comparing the C 18 Q(2- 1) data of IBADBO2I with 
the 1.3mm continuum map of MAOll we measure a reduction 
of the C ls O abundance (averaged along the line of sight) by 
a factor of ~ 3.5 toward the center of the envelope, assuming 
an optically thin C ls O line 1 . There is thus no evidence for any 
C I8 (and thus CO) desorption from the dust grains in the inner 
envelope at the 30m resolution (800 AU in radius). This indi- 
cates that potential CO desorption due to shocks produced by 
the outflow does not dominate in the central beam. Although 
the overall morphology seen in Figs. [2] and |4^ suggests that 
the N2FF abundance may reach a minimum along the outflow 
path, the outflow is sufficiently well collimated (cf. Fig. 1 of 
AMB99) that its effect should be negligible in the equatorial 
plane of the envelope. 

1 Using our C 17 (l-0), C 18 (l-0) and C 18 (2-l) data, we estimate the 
C 18 (2-l) optical depth to be 0.7-1 toward IRAM 04191. 



We suggest two other possible explanations for the disap- 
pearance of N 2 H + in the IRAM 04191 envelope: 

1) N2 may condense onto the grain surfaces. This could be 
either due to a higher N2 binding energy on polar ice man- 
tles than the value of 750 K usually a dopted (see Model C of 
lAikawa et alJ200lllBergin et alJ2 002). or to an enhanced grain 
chemistry transforming N2 into less volatile species (Y. Aikawa 
2003, private com municat i on). S uch a strong N2 depletion was 
suggested bv lCaselli et alJ ll2003h to account for the high H2D + 
abundance measured in LI 544. 

2) Deuterium fractionation may enhance N2D + over N2FU. A 
high-degree of deuterium fractionation has indeed been ob- 
serv ed in dense co r es (e.g . |Turnerl200lHBacmann et all2003l) . 
and iRoberts et ail d2003l) predict that the N 2 D + /N 2 H + abun- 
dance ratio should even exceed 1 at densities of ~ 10 6 cm -3 
when they include HDj and in their chemical model. N2D + 
may thus be a better tracer of dense gas than N2I-U. 

The disappearance of N2FU from the gas phase suggests 
that all species bearing heavy ele ments also vanish at high den- 
sity in prestellar cores. However, IWalmslev et alJ (2004 find 
that, under such conditions, the ionization degree of the gas is 
not significantly changed owing to the enhanced abundances of 
H + , H3 , and . Moreover, they argue that the ambipolar dif- 
fusion timescale is dominated by the coupling of the neutrals 
to the charged grains at these densities. The dynamics of proto- 
stellar collapse may thus be essentially unaffected by the strong 
molecular depletion in prestellar cores. 

Acknowledgements. We would like to thank Y. Aikawa, M. Gerin, G. 
Pineau des Forets and M. Walmsley for enlightening discussions about 
the chemistry. We are grateful to F. Gueth, J. Pety and R. Neri (IRAM) 
for their help with the Plateau de Bure observations. We would also 
like to thank F. Motte for her advice about the continuum data. 

References 

Aikawa, Y., Ohashi, N., Inutsuka, S.I., Herbst, E., & Takakuwa, 

S. 2001, ApJ, 552, 639 
Aikawa, Y, Ohashi, N., & Herbst, E. 2003, ApJ, 593, 906 
A ndre, P., M otte, F, & Bacmann, A. 1999, ApJ, 513, L57 

ilAMB99h 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, in 

Protostars and Planets IV, ed. V. Mannings, A.P Boss, & S.S. 

Russell (Univ. of Arizona Press, Tucson), p. 59 
Bacmann, A., Lefloch, B., Ceccarelli, C, Castets, A., 

Steinacker, J., & Loinard, L. 2002, A&A, 389, L6 
Bacmann, A., Lefloch, B., Ceccarelli, C, Steinacker, J., 

Castets, A., & Loinard, L. 2003, ApJ, 585, L55 
Belloche, A Andre, P ., Despois, D., & Blinder, S. 2002, A&A, 

393, 927 JBADB02I) 
Bergin, E.A., Alves, J., Huard, T, & Lada, C.J. 2002, ApJ, 570, 

L101 

Caselli, P., van der Tak, F.F.S., Ceccarelli, C, & Bacmann, A. 

2003, A&A, 403, L37 
Doty, S.D., & Neufeld, D.A. 1997, ApJ, 489, 122 
Greve, A., Kramer, C, & Wild, W. 1998, A&AS, 133, 271 
Guilloteau, S., 2001, in Proceedings from the IRAM Millimeter 

Summer School 2, ed. A. Dutrey, p. 219 



A. Belloche and P. Andre: Disappearance of N2H + from the Gas Phase in the Class Protostar IRAM 04191 



5 



Motte, R, & Andre, P. 2001, A&A, 365, 440 (MA01) 
Motte, R, Andre, P., & Neri, R. 1998, A&A, 336, 150 
Roberts, H., Herbst, R., & Millar, T.J. 2003, ApJ, 591, L41 
Shematovich, V.I., Wiebe, D.S., Shustov, B.M., & Li, Z.Y. 

2003, ApJ, 588, 894 
Tafalla, M., Myers, PC, Caselli, P., Walmsley, CM., & 

Comito, C. 2002, ApJ, 569, 815 
Turner, B.E. 2001, ApJS, 136, 579 

Walmsley, CM., Flowe r, D.R., & Pi neau des Forets, G. 2004, 

A&A, in press |astro-ph/0402493| 
Wiesemeyer, H. 2001, in Proceedings from the IRAM 

Millimeter Summer School 2, ed. A. Dutrey, p. 71 



